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Abstract 

We derive the fully backreacted bulk solution dual to a boundary superfluid 
with finite supercurrent density in AdS/CFT. The non-linear boundary 
hydrodynamical description of this solution is shown to be governed by a 
relativistic version of the Tisza-Landau two-fiuid model to non-dissipative 
order. As previously noted, the phase transition can be both first order 
and second order, but in the strongly-backreacted regime at low charge q 
we find that the transition remains second order for all allowed fractions of 
superfiuid density. 



1 Introduction 



A holographic superfluid PQ[2] is a charged, asymptotically AdS black hole, which below 
a certain temperature Tc spontaneously develops a scalar condensate. The relevant 
instability of charged scalars in AdS-black-hole backgrounds was first pointed out in |3]. 
In many contexts one can also use such systems to explore the physics of holographic 
superconductivity. In this work we focus on the superfluid aspects by deriving the 
boundary superfluid hydrodynamical equations that govern them. 

According to the AdS/CFT correspondence [1], asymptotically AdS black holes are 
dual to thermal states of a large- boundary conformal field theory If we re- 
strict ourselves to considering low-frequency fluctuations of this field theory with large 
wavelengths, we expect an effective hydrodynamic description to be valid [6j. Said 
differently, the hydrodynamical description is applicable for configurations of the field 
theory which are locally in thermodynamic equilibrium and which vary from local 
patch to local patch such that all gradients remain small. It was shown in [7] that 
such a field-theory gradient expansion can be implemented in AdS/CFT as a pertur- 
bative scheme of solving Einstein's equations in the universal subsector of gravity with 
a cosmological constant. Charged branes, and thus chemical potentials for conserved 
global currents in the boundary theory, were introduced by [8l |9]. The setup was sub- 
sequently generalised in [lOj to include effects a bulk U{1) field inducing boundary 
magnetic fields. 

A qualitatively rather different picture should arise for systems that allow for the 
breaking of a global U{1) symmetry thus giving rise to superfluidity. The hydrodynamic 
limit of such a system is given by the Tisza-Landau pTl [T2] two-fluid model, which 
describes the superfluid as a two-component mixture of a normal fluid with density p„ 
and velocity and a superfluid component with density ps and flow velocity v^. 

In this paper we find solutions of Ad gravitjj^ which we argue to be dual to holographic 
superfluids with finite supercurrent density and we show that their hydrodynamical 
description is captured precisely by a relativistic version of the Tisza-Landau two- 
fluid model. We do this by constructing asymptotically AdS solutions that break 
the bulk U{1) symmetry with a charged condensate at finite chemical potential. The 
use of charged AdS black holes to study field theory at finite chemical potential was 
considered early on in albeit in the context of global AdS. Our solutions are gauge- 

""^The spacetime coordinates {x^^} split into {r, x^}, where r is the holographic direction and we 
refer to {x^^} as the field-theory or boundary coordinates. 
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equivalent to situations where the phase of the bulk scalar, given the interpretation 
of the field-theory Goldstone mode, has non-trivial dependence on the field theory 
directions. To non-dissipative order in the fluid expansion this dependence is linear so 
that the Goldstone mode has a constant gradient. 

Some properties of holographic superfluids, both with s-wave and p-wave conden- 
sates, such as its various sound modes etc., were explored in [HI [15], HEl [TTl [181 [iHl [201 
m\ . Our treatment is complementary to theirs in two ways: Firstly we improve on the 
state of the art by including the backreaction consistently, that is we solve the full set 
of coupled Einstein-Maxwell-scalar equations. Secondly, our approach, akin to that of 
[7] for the normal fluid case, derives the non-linear boundary fluid dynamics, including 
the 'Josephson equation' from bulk gravity, rather than taking them as a given and 
then fitting transport coefficients to the postulated model. 

The paper is organised as follows. Immediately following this paragraph we give a 
brief overview of the analytical results of this paper. In section 2 we review and extend 
the by now well-known holographic superconductor. In section 3 we describe the gravity 
setup corresponding to introducing a non-zero supercurrent on the boundary, perform 
its holographic renormalisation and define the appropriate statistical grand partition 
function. In this section we also derive the boundary stress tensor and U{1) current. 
Section 4 contains our numerical method and results and section 5 gives a discussion 
of the results and some open issues. Appendix A contains the bulk equations in their 
full glory. 

1.1 Summary of results 

Here we summarise the findings of the rest of the paper. We emphasise that the equa- 
tions of fluid dynamics, as well as the form of the constitutive relations are exact and 
derived analytically. The hydrodynamical equations following from our bulk analysis 
to non-dissipative order can be summarised as follows 

(e + P) u^Uu + Pr]^^ + ixpsV^Vu 

J PL = PnU^ + PsV/, , 

d^J^ = 

d^T^" = (1.1) 

where the total charge density is p = pn + Ps and we have set a possible external applied 
field that can appear on the right-hand side of the stress-tensor conservation equation. 
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to zero. The quantity has the interpretation of normal fluid velocity, whereas is 
the superfluid component. These satisfy the constraint 

= -1 (1.2) 

also known as a Josephson equation. The rest of this paper explains the bulk origin of 
these equations and numerically constructs the dual solutions. We also use our non- 
linear solutions to study certain physical properties, such as the order of the phase 
transition, numerically. 



2 Isotropic holographic superfluids 

In this section we review and (slightly) extend the work of [H E] on the equilibrium 
properties of holographic superfluids. This section also serves to introduce our notation 
and conventions: rather than using standard Schwarzschild-like coordinates, we use 
ingoing Eddington-Finkelstein coordinates and introduce a boost along the direction 
of the (flat) horizon. The class of solutions we are constructing in this section are thus 
generalisations of boosted black branes allowing for a non-trivial proflle of a complex 
scalar fleld. 

A simple action that is sufficiently general for our aims is given by 



^ -FmnF''^ -vm- w? 



d^x, (2.1) 



supplemented by the Gibbons-Hawking- York boundary term. Here t/^ is a complex 
scalar with charge q under the U{1) gauge fleld 

Dil) = di)- iqAijj (2.2) 
and the potential V^dV'l) is chosen to consist of just a mass term 

Vm) = -h^\\ (2.3) 



where £ is the AdS length. Different choices of the potential j22t 123] are possible. The 
work of [16] used flrst-order perturbation theory to study the effects of a {ipl'^ term 
on the various superfluid sound modes. Embedding in M-theory [211 EH] or string 
theory |2Bj dictates a speciflc form of the scalar potential and it would be interesting 
to investigate their hydrodynamics. 
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As we are interested in black holes with flat horizon sections, we propose an ansatz 
of the form 



^"^static — —"^hundx^dr — fu^u^dx^dx'^ + ^Afj_,ydx^dx'^ 

f(l)h 

A = -(pu^dx^ 

ijj = l^le^^" := ^e^«". (2.4) 



We have chosen to write everything covariantly in the boundary directions {x^} by 
introducing the boost velocity which satisfies 7^m^u'^ = — 1- The quantity 7^ : = 
limj._j>oo is the metric on the boundary 



dsi,=^^JxV. (2.5) 
Finally, we have the spatial projector 

= 7~. + Uf,u, (2.6) 

At this point this is nothing more than a boosted version (by u^) of the holographic 
superconductor solution of [2J written in Eddington-Finkelstein coordinates. 

If we rescale A^ ^ ^ and ip ^ the rescaled Maxwell-scalar sector is down by 
a factor of q^^ with respect to the Einstein-Hilbert action. Taking the limit g — )■ 00 
thus eliminates any backreaction on the metric. Gravity backreaction is physically 
important not only for low charge q but also at low temperatures. In this paper we are 
able to work away from the probe limit. For completeness we demonstrate that our 
results do indeed go over into the various perturbative results known in the literature 
in the appropriate limits. 

2.1 Exact black brane solution 

The normal phase of the superfiuid which has vanishing order parameter on the bound- 
ary is given by the standard charged RNAdS black brane 

f = h 1 — — 

4> = hoc ifi ^j' h = hoo, V = 0. (2.7) 

The metric on the boundary is now 

dsl^ = -hl^dv^ + Sijdx'dx^ . (2.8) 
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In order for A to be a well-defined one form at the horizon (r = r_|_), we must set 

P=^. (2.9) 

This can be seen, e.g. by going to Kruskal coordinates. In this gauge fi can be identified 
as the chemical potential of the boundary theory. Also p is the charge density of the 
system 

*F = pvol2 (2.10) 



with V0I2 the divergent two-volume of the fiat spatial sections of the geometry. The 
Euclidean solution can be thought of as a thermodynamic ensemble at temperature 

We will soon generalise the metric ansatz, but the expression for temperature will 
remain unaltered. 

2.2 Holographic stress-energy and current 

It is well known that the naively defined boundary stress tensor in asymptotically AdS 
space needs to be renormalised by introducing suitable counterterms [27]. For our 



system (2.1), the desired stress tensor and current follow from varying the action 



Sren = ^bulk + ^ j + ^ + j [ ) £x , (2.12) 



SttG 



where 7^1, is the induced metric on the boundary. In addition to the standard gravita- 
tional terms there is an additional counterterm for the scalar field. This choice is valid 
provided that one fixes the coefficient of the leading behaviour of ip to zero. We 
restrict to this boundary condition from now on. Alternatively, we could have fixed the 
coefficient of the behaviour by adding a different counterterm for the scalar field tp. 
The mass of the field ip places it in the range where both behaviours can be thought of 
as vevs. For a more extensive discussion, involving both possible choices of boundary 
conditions, see e.g. [21 ES]- order to evaluate the stress tensor and current, we need 
to determine the asymptotic expansions of the fields. The near-boundary expansion of 
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the equations in appendix \K\ yields 

h = hn 



1 ^ + 



/ = 


, 2 STiGfe 


= 


hoo\l^ 

y r 


i = 





+ 
+ ■ 



(2.13) 

Furthermore, we can set h^o to unity by invoking a scahng symmetry and we will do 
so wherever convenient. However, it is important to keep in mind where it occurs 
formally when we compute boundary variations later on and variations of h^o translate 
into variations of temperature. 

It now follows that the current and stress tensor 



Jn 



2 55'ren 
1 (^iSron 



SvrG 



are finite in the limit r — )■ oo. An explicit calculation reveals that the background 
solution reproduces the stress tensor and current of an ideal conforma|^ relativistic 
fluid, which we can see from the leading order behaviour 



r 

lim -T, 



r 

hm -J„ 



We can also determine the total energy of the solution, 

/v 
d^x-Too = evoh, 



(2.15) 



(2.16) 



where the total (divergent) two-volume V0I2 comes from the integral over spatial direc- 
tions. 



2.3 Thermodynamics 

Having already determined the counterterms, we can use standard Euclidean techniques 
[28] to find the thermodynamics of this system. We want to calculate the Euclidean 

^One can show that TrT^^ cx ^1^2, so that the fluid satisfies the tracelessness requirement of 
conformal invariance if and only if either ^1 or ^2 is zero. In the dual field theory this corresponds 
to not deforming by the relevant operator O^. Of course finite temperature and/or finite chemical 
potential then break conformal invariance. 
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action Ie, which takes on the role of the Gibbs free energy for the holographic superfiuid 
in the saddle point approximation. The results in this section are standard so we will 
give few details on their derivation. We take 

V = -IT , Ie = -iSren ■ (2-17) 

and use the fact that the gravitational Euclidean partition function defines a thermo- 
dynamic potential via 

\nZ = -l3n. (2.18) 

Following the same kind of manipulations as in [251 El] we find that the manifestly 
finite Euclidean action gives us the thermodynamic potential 

Ie = /3fi(/i, T) = ^ (e - /ip - Ts) , (2.19) 

where s is entropy density (the entropy S* = ^ is defined as usual in terms of the 
horizon area). From the on-shell variation of the action we learn that the action is 
stationary for fixed temperature and fixed chemical potential fi. Furthermore, 
stationarity of the action requires that either = or ^2 = 0. Thus we have the first 
law 

= T5s + fi5p . (2.20) 

We will now generalise this system by adding an additional chemical potential, related 
to the Goldstone field of the broken global U{1) of the boundary theory. This has the 
physical interpretation of allowing for a finite supercurrent to flow in the broken phase. 
The bulk manifestation of this is a non-zero value for the spatial components of the 
vector field, best quantified by saying that the gauge invariant expression — d^^a 
develops a non-trivial value in the bulk. 



3 Introducing a Supercurrent 

Let be a constant vector of unit magnitude, such that 

7^X«- = (3.1) 
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This means that in a frame where = 1 , u = 0, we can write ra" = and the spatial 
components are given by n, a unit two vector. We consider the bulk ansatz 

^2 2Ch 

ds^ = (iSgta^tij, + — (2 Cu(^n„) - Bn^nJ) dx'^dx" + —^n^dx'^dr 

A = {-(t)u^ + $72^) dx^" - ^dr 

%l) = ^e*"". (3.2) 
The additional functions in this ansatz depend only on the radial direction. The con- 



dition (3.1) can be imposed without loss of generality in the sense that any constant 



space-like vector defines a bulk metric that is diffeomorphically equivalent to (3.2). 



So, by splitting up such a space-like vector into components parallel to u'^ and compo- 



nents normal to we can arrange for (3.1) to hold. 



For specific calculations we will always work in a gauge where the scalar ip is real. 
However, notice that such solutions are gauge-equivalent to ones where the phase a 
has a boundary dependence that is linear in boundary coordinates. To see this, one 
applies a gauge transformation 

A^^ A^ + d^ix^X^), a^a + x^K. (3.3) 

Let us define the combination 

«M = -4M-VMa- (3.4) 

Then the quantity 

= lim {A^ - d^a) (3.5) 

is gauge invariant with respect to the bulk U(l) and therefore an operationally conve- 
nient measure of the magnitude of the boundary Goldstone mode. We note that 

a^a^^ = -fi' + C (3.6) 

gives us a gauge-invariant definition of the boundary supercurrent density. The lat- 
ter part of this paper will be concerned with numerically solving the Einstein-Scalar- 
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Maxwell equations 



Vm(^VO = (3.7) 

with the above ansatz. They are given explicitly in appendix A. However, even without 
numerics, we can learn a lot about the system by analysing its behaviour in an analytic 
near-boundary expansion. Let us turn to this analysis now. 

3.1 Asymptotic analysis 

The equations allow the following asymptotic expansion (the remaining fields have the 
same form as in ( 2.13[ ) ): 



^3 3?^"^ 



/ = ft--— 3 + 



In order to keep the asymptotic boundary metric in Minkowski form, we demand that 
the non-normalisable parts of B and C vanish. Thus we restrict attention to solutions 
with 

g(o) ^ (3(0) ^ _ (3 9) 



3.2 Free energy and First Law 

We can now identify the relevant thermodynamic potentials and variables of the bound- 
ary theory from our asymptotic bulk analysis. In order to relate the Euclidean action 
to a thermodynamic partition function we first have to write the bulk part of the action 
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as a total derivative. This can be done in two different ways: 



5, 



bulk 



(11) 
bulk 



167rG 

2/3vol2 



oo /2 

r 



-99 



oo 



dr 



(3.10) 



The first expression gets contributions only from the upper boundary of the integral, 
whereas the second expression receives both horizon and boundary contributions. 

Let us consider the on-shell variation of the Euclidean action. This receives no 
contribution from the lower boundary of the integral. From the upper boundary we 
find for variations that preserve the (conformally) Minkowski metric on the boundary: 



For fixed /3, we have Sh^o = - 



OS 
ST 



Pvoh [{e - fip) 5hoo- pSfi + JsSfis 



(3.11) 



, ^, so that the action is stationary for fixed li,Us,T. 
Putting all this together we learn that the Euclidean action defines a grand canonical 
ensemble in terms of temperature T and chemical potentials p and ps- 
Using 



rbulk 



counter 



we find the quantum statistical relation 
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i3(3) 



where P is pressure, following from the first form of the action, and 

uj'^^^\fi,fis,T) = 2P -Ts - fip + fisJs 



(3.12) 



(3.13) 



(3.14) 



from the second form. Equating the two different expressions for the thermodynamic 
free energy gives rise to a Smarr-Gibbs-Duhem relation: 



3P = Ts + pp — psJs ■ 



Using (3.14) and (3.11) we can deduce the first law in the form 

e-2P 



6uj = — \ s + 



T 



psJ s 



6T - p6p + JgSps 



(3.15) 



(3.16) 



Note that via (3.6) and (3.8) the third term involves a variation of the norm of the 



boundary Goldstone field. Using the relation u = —P we can rewrite this in the form 



6P = s6T + pJp-^6{{a^f) 



(3.17) 
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where we defined the normal and superfiuid density as 



P = Pn + Ps, with p^ = — (3.18) 

Ps 

The conjugate function of temperature, 

s = s + , (3.19) 

is the thermodynamic entropy, which gets contributions not only from the horizon 
entropy, but also from the hair of the black hole. We will see in the next section that 
these contributions in fact cancel each other, so that the final expression of the entropy 
density is that of the black-hole horizon. Thus the pressure P (T, fi, (a^)^) is the same 
thermodynamic potential as that appearing in the context of relativistic field theory 
in [291 EQ]- Here we gave a derivation from bulk gravity, which via the AdS/CFT 
correspondence is dual to a strongly-coupled relativistic field theory on the boundary. 
We can rewrite the Smarr-Gibbs-Duhem relation as 

P + 6 = Ts + fxpn. (3.20) 

It is obvious from this definition that e is the Legendre transform of pressure with 



respect to temperature and chemical potential. From (3.15) and (3.18), it is related to 



the zero-zero component of the energy-momentum tensor via 

e = e-pps. (3.21) 

3.3 The energy momentum tensor and current 



Using the renormalised expressions (2.14) we can deduce the form of the boundary 



stress tensor and current following from the bulk ansatz (3.2). We find 

T^u = {e + P)Uf,u^ + Pr]f,^-B^'^^nf,n^ + 2C^^^u^^n^) 
JfM = pu^ - JsTif, . (3.22) 

Notice that the stress tensor is traceless on account of the orthogonality relation be- 
tween and and the definition of pressure in terms of e and B^^\ One might be 
a little surprised by the occurrence of a cross term between and n^. To elucidate 
this point, let us compare our form of the constitutive relations with those of [HI [29] . 
These in turn were shown to be equivalent to the original formulation of Israel, Carter, 
Khalatnikov, and Lebedev (see references in [JSUSH])- 
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Let us define a new superfiuid velocity via 



n„ 



g(3) 



with 



(3.23) 



In the new frame we get a stress tensor and current 

p(3)]2 



e + 



J„ 



^3(3) 

C(3) 



+ P u^u^ + Pr/, 



J. 



■u^ + J, 



C(3) 



(3.24) 



The coefficient of is the superffuid density ps- From section 3.2 the entropy density 



in the ffuid rest frame is s. To get to the arbitrary reference frame of (3.24) we must 
apply the boost u^. From this, it follows that the current su^ describes the flow of 
entropy. As in [llj, one can constrain the form of the constitutive relations by requiring 
a vanishing divergence of the entropy current 



This requires that 



with 



ps 



J. 



C(3) 

m 



(3.25) 



(3.26) 



Finally, we should compare this expression for the superfluid density to that in Eq. 

C(3) 



(3.18). Evidently we must have 



p_ 

Ps 



(3.27) 



Interestingly, using (3.27) and (3.26) in (3.19) we see that the additional contribu- 
tions to the entropy cancel against each other, so that s = s. We just derived these 
constraints among the boundary data from a careful analysis of the thermodynamics 
following from a semiclassical treatment of the bulk quantum partition function. The 



conditions (3.26) and (3.27) have to be satisfied if the boundary is to have a sensi- 



ble thermodynamic description and indeed they imply that the entropy s is in fact 
the thermodynamic entropy s of the black hole. However, we cannot simply impose 
these constraints on the asymptotic data of any generic solution. There is not enough 
freedom in the equations to do this. We therefore need to check that the full bulk so- 
lutions, once appropriate regularity conditions have been given at the horizon, satisfy 
these constraints automatically. We will see below that both relations between asymp- 
totic data coefficients are satisfied by our numerical solutions to very high accuracy 
(see figure |4]). We should see this result as a non-trivial consequence of the interplay 
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between horizon regularity and boundary fluid dynamics that can only be seen at the 
full non-linear level. 

The description of superfiuids we have derived in the preceding sections is identical 
to that proposed by Landau pJj, and later generalised to the relativistic domain by 
Khalatnikov et al., |31j . 

Note that when going on to include dissipative terms, i. e. terms of higher order in the 
derivative expansion, the conservation of the entropy current should be replaced by the 
requirement that its divergence be non-negative. A positive entropy current has been 
constructed for general fluid-dynamical solutions without bulk symmetry breaking in 
[32j and it would be interesting to see if such a construction could also be carried out 
in our case. The analysis in this paper establishes semi-analytically - in the sense that 
the horizon-to-boundary map had to be constructed numerically - that such a current 
exists at non-dissipative order. 

Finally note that in terms of v^, the gauge invariant quantity evaluates to 

a';: = -fiv,. (3.28) 

By choosing a different gauge in the bulk, we can make direct contact with the standard 
treatment, given for example in [21] • Thus, if we take d^a = /if^, the Josephson 
condition becomes 

fi + M^(9^a = , (3.29) 
where now /i is the chemical potential appearing in [29]. Note that we had to use 



the constraints (3.26) and (3.27), which we have confirmed numerically, in order to 
connect the quantities ra^ and consequently to the superfiuid velocity. Analytically 
these relations follow from requiring a positive semi-definite divergence of the entropy 
current. However this property of the divergence should follow from the second law 
of black-hole thermodynamics and establishing this law in the present context would 



constitute an analytical proof of the identifications of (3.26) and (3.27). 



3.4 Local bulk symmetries — > Boundary hydrodynamics 

One use for the solutions we have constructed is as a homogeneous starting point 
for a hydrodynamic gradient expansion of the form pioneered in [7]. Although many 
interesting physical properties have been obtained in (ordinary) perturbative studies 
[m Uni [in], these perturbative solutions do not serve as a good starting point for a 
gradient expansion and we comment on this further in the discussion section. We hope 
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to return to the dissipative order of the gradient expansion in the future. However, 
even without exphcitly carrying out the higher orders of this expansion, we can already 
deduce the hydrodynamic relations these will satisfy from general symmetry principles. 

Among the Einstein equations there are constraint equations, which are first order 
in radial derivatives, and dynamical equations, which are of second order. In the 
context of a gradient expansion of the universal gravity subsector of AdS/CFT it 
was shown in f7] that these constraint equations at order n imply that the boundary 
stress tensor and U{1) current satisfy the equations of fluid dynamics up to order n — 
1. However, these equations in fact follow from elementary symmetry considerations. 
Namely, boundary current conservation follows from the U{1) gauge symmetry in the 
bulk, and the remaining equations, involving the stress-energy tensor, can be derived 
by considering an infinitesimal diffeomorphism generated by a vector field lying parallel 
to the boundary directions. 

The resulting equations take the form of a conservation equation for the current and 
an equation of continuity for the stress-energy tensor: 

9^ = , S^T'^^ = ^xJ^ (3.30) 

where T is the applied boundary field strength tensor (to be distinguished from the 
bulk field strength F). The bulk gradient expansion then serves to derive the form 
of the conserved stress and current order by order. In this work we have determined 
all non-dissipative contributions. It is a straightforward extension to go on to include 
dissipative terms. 

From the bulk analysis we have also learned that 

v% = -l, (3.31) 

which gives us the 'Josephson equation' governing the dynamics of the Goldstone mode. 
From a field-theory perspective it has been elaborated upon in p9]. Together with 



the constitutive relations (3.22), the equations (3.30) and (3.31) form the full set of 
equations of conformal relativistic superfiuid hydrodynamics. This is as far as the 
analytical results can take us. We will now turn to the numerical part of this paper. 

4 Numerical Results 

We wish to solve the differential equations of appendix |A] subject to the following two 
requirements 



14 




Figure 1: (colour online) The condensate {O^) as a function of reduced tem- 
perature. For larger values of the charge q the condensate becomes multi-valued 
beyond a critical fraction of superfluid Cc indicating that the phase transition 
is now first order. The normalisation of the condensate is chosen for ease of 
comparison with the probe calculations of 114] - 



1. There is a regular horizon in the bulk at r = r+. 

2. The asymptotic boundary metric remains conformally flat. 

This is the minimal set of physical requirements that ensures that the boundary theory 
is dual to a locally equilibrated fluid living on Minkowski space. From the second 



requirement we already deduced the condition (3.9) on the asymptotic expansion (3.8). 
It is convenient to choose units (that is using a scaling symmetry) to fix r+ = 1. 
While the second condition has already been incorporated in the asymptotic boundary 
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expansion, the first requirement fixes the near-horizon expansion to take the form 



f{r 
h{r 

C[r 
0(r 
$(r 



U{r - r+) + ■ • • 
/i+ + ■ • • 

C+(r-r+) + --- 
0+(r -r+) 
$+ + ■■■ 



with 



where C-( 



/+ = hi {3 + WnG^l) - 47rGf(P+ , (4.1) 
$+, ^+} are the independent horizon data and /+ and all higher- 



order coefficients are uniquely determined by them. Regularity demands that C van- 
ishes at the horizon in such a way that j remains finite, as can be seen for instance by 



inspecting the n^dx^dr term of the metric in (3.2). Similarly 0(r+) = for regularity. 
In practise we determine the expansion to a high order and use it to set boundary 
conditions near the horizon. 

The total differential order of the system of ODEs we are solving is twelve. After 
taking account of all scaling symmetries we are left with eight more data at infinity, 
, /is, Jsj ^oo;£^}; giving a total of 14 independent pieces of data. We will 
thus get a two-parameter family of solutions. Of course we are free to choose any two 
non-confiicting parameters, but a natural and convenient choice is the temperature T 
and the superfiuid chemical potential /i^. Because of the underlying conformal symme- 
try, all physical quantities will depend on the dimensionless ratios ^ and C = We 
will refer to t = 



T 



qfj. 



as the reduced temperature and ( as the superfiuid fraction. 



4.1 The condensate 

The asymptotic behaviour of the field if) sets the expectation value of the symmetry- 
breaking dimension-two operator in the boundary theory. 

{O^) = 726 (4.2) 

This operator tends to condense at low temperatures. By numerically constructing 6 
we determined the graphs shown in figure [T] illustrating the dependence of the con- 
densate on reduced temperature and superfiuid fraction. There is always a critical 
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0.00 0.01 0.02 0.03 0.04 0.05 0.00 0.01 0.02 0.03 0.04 0.05 

Figure 2: (colour online) The quantity ps/p as a function of reduced tem- 
perature. This fraction approaches unity for q > 1, but for q = 1 there always 
remains a normal component, even at very low temperatures. Some of the curves 
are multi-valued as a consequence of the first-order phase transition. In those 
cases only the upper branch below the critical temperature tc{C) is physical. 



superfluid fraction above which the system is forced into its normal state. The same 
behaviour was found in probe-hmit calculations [Ml [I5], as well as for rotating su- 
perfluids in [33] • Furthermore, for sufficiently high charge q of the field dual to the 
condensing operator there is also a special point on the phase diagram below which the 
superfiuid transition is first order and above which it is second order. This behaviour 
was already found in the probe limit calculations. However, we find that for lower 
values of q the transition never becomes first order. This is clearly a result of strong 
gravitational backreaction. In order to conclusively show that the transition is first 
order or second order we also computed the free energy as a function of temperature 
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for fixed (. 



4.2 Free energy and superfluid density 

Figure [3] show numerical graphs of the thermodynamic free energy 

(4.3) 



/i^vol2 dT 



at fixed superfluid fraction. We show curves for q = 1,2. For sufficiently high q and 
larger fractions (, we clearly see that the smooth second-order behaviour goes over into 
the swallow-tail cusp behaviour characteristic of first-order phase transitions, but for 
q = 1 the free energy remains smooth and the transition second order up until the 
critical value (c above which no superfluid phase exists. The free energy higher values 
of q looks very similar to the q = 2 case and we do not show the corresponding plots 
here. The bulk charge q appears in the boundary theory as a normalisation factor of a 
ratio of three-point and two-point functions [3J, so from this perspective there appears 
to be no need to restrict its value to be above a certain treshold or indeed to be integer 
valued. We use integer values in our numerics only for convenience, as is the case in 
most of the pertinent literature. Figure |2] presents plots of the behaviour of p„ and ps as 




Figure 3: (colour online) The free energy at fixed C, as a function of reduced 
temperature for q = 1,2. The unbroken branch is shown as a dashed line. For 
the lowest q the broken branch joins the unbroken branch smoothly for all allowed 
values of C,. For q = 2, we see the swallow-tail behaviour characteristic of first- 
order transitions for fractions above a critical value Qc- The behaviour for q > 2 
is very similar and we do not reproduce these plots here. 
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Figure 4: (colour online) Figure demonstrating to numerical accuracy (in this 
case the numerical error is of the order of 10~'^) the constraints among asymp- 
totic data following from horizon regularity. Both constraints are illustrated based 
on the data points we used to generate the plots at q = 8 and C = i • 



a function of temperature. The double-valued nature of the graphs is a consequence of 
the first-order phase transition. At Tc the system spontaneously develops a finite value 
of the condensate and immediately jumps to the thermodynamically preferred branch. 
Again we emphasise that the low-temperature behaviour involves strong backreaction 
and can therefore not be determined in the probe limit. 

Figure |4] demonstrates the validity, to numerical accuracy, of the relations between 
the two expressions for the superfluid density, Eqs. (3.26) and (3.27), that we derived 



in section 3.3 from basic thermodynamic considerations. 

Finally, in figure [s] we show that the critical temperatures of [H] are recovered for 
large q for appropriately rescaled quantities. Note that even though q = 1,2 are far from 
the probe limit, q = 8 already gives a close approximation to the critical temperature 
in the probe limit. 

We shall conclude with a brief discussion of the relevance of our approach to de- 
termining higher-order transport properties of strongly coupled relativistic superfluids 
and mention some open questions. 



5 Discussion 



If one attempts to formulate a gradient expansion of a hydrodynamical system all 
transport coefficients at non-zero derivative order are subject to an ambiguity stem- 
ming from the freedom to shift zeroth-order quantities by gradients. In the gravity 
setup these shifts manifest themselves as homogeneous modes of the bulk equations. 
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Figure 5: (colour online) Temperature, at which the marginal mode occurs as a 
function of superfiuid fraction (. The curves for ascending q = 1,2, 8, 16, 32, 100 
are shown from bottom to top. The vertical dashed lines indicate fractions of the 
condensate of C, = 1/4,1/3,2/5,1/2 and the horizontal dashed lines are visual 
aids in order to identify the critical temperatures for the q = 100 curve. 

One can resolve this ambiguity by fixing appropriate frames, a popular choice being 
the Landau frame. Alternatively one can concentrate on a set of scheme- independent 
quantities. However, both approaches show that the set of scheme-independent quan- 
tities (transport coefficients) at first order, and therefore in turn all higher orders, 
depends on the number of independent thermodynamic quantities present already at 
zeroth order. This means that, for a well-defined gradient expansion that includes 
non-trivial transport coefficients associated to the superfiuid component, we must in- 
corporate the superfiuid velocity in the homogeneous zeroth-order solution. We have 
constructed such a supercurrent solution in this paper. Incorporating the supercurrent 
only perturbatively at zeroth order p5] means that the hydrodynamical expansion 
must be a double expansion, both in slow gradients and small superfiuid density. This 
is appropriate for determining physical quantities, such as various sound speeds away 
from low temperature or large condensates [171 HE] , but we would expect the method 
to break down at some point, when the superfiuid component strongly backreacts on 
the metric. It is thus an interesting project to determine systematically, using our 
techniques, the new transport coefficients associated with the superfiuid component to 
first and perhaps second order. 

This work assumes that classical (super) gravity is a good approximation to the 
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Figure 6: (colour online) Curvature at the horizon for different fractions of 
superfiuid. The solid lines depict the Ricci scalar as a function of temperature 
and the dashed lines the Kretschmann scalar Kr = ^ RmnpqR^^^^^ : whose 
normalisation N is chosen so that Kr = N at tc{C = 0) . Both quantities diverge 
at very low temperatures, indicating that the supergravity approximation breaks 
down. 



underlying string theory. This is true so long as the geometric curvature scales are 
small and thus the dual gauge theory is strongly coupled. Figure |6] demonstrates that 
the assumption that curvatures remain small is not valid at low temperatures close to 
the horizon for g = 1, 2 (again, the qualitative behaviour for g > 2 is similar to g = 2 
and thus not shown). This was demonstrated previously for zero superfiuid fraction in 
[211 |25l 131] . It was also shown that adding higher-order terms to the scalar potential, 
as dictated by M-theory, can cure this problem. We expect the same to be the case 
here. 

Furthermore we demonstrated that the holographic superfiuid admits an entropy 
current, associated with its normal component, with positive definite divergence. Fun- 
damentally, in the AdS/CFT setup, this of course follows from the second law of black 
hole thermodynamics and thus there should be a general proof to all orders in the gra- 
dient expansion. Such an an entropy current was developed in full generality for the 
gravity subsector in [32] and it would be interesting to see if a similar analysis can be 
made for our solutions. We have reason to be optimistic about this, since the entropy 
of the field theory at zeroth order is simply given by the classical entropy associated 
with the black-hole horizon. 

Finally, the analysis of this paper focuses on bulk solutions of Ad gravity, corre- 
sponding to thin films of superfiuids. It would be interesting to repeat the numerical 
analysis in one higher dimension corresponding to relativistic 3 + 1 dimensional su- 
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perfluids which are also relevant for certain parts of the QCD phase diagram]^ Many 
of our analytical results are valid in any dimension and the remaining ones can be 
straightforwardly extended. 
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A The bulk equations 

Here we reproduce the bulk equations governing the holographic superfluid. The equa- 
tions for /(r) and h{r) are first-order ODEs, the rest are second order. 

/ equation: 

-r'^f [rB' + 46-4) /' = 

Ar^fB' + Ar'^h^^^ {-Bf + + /) - 12r^ (/ - h"^) {-Bf + + f) 
+2rVC'^ {-Bf + C^ + f)+ r^(B - l)f4>'^ - 2q^r^{B - l)h^^^(j)^ 

+4q^r^Ch^^^(f>'P - 2q^ f h"^ + r'^f<^>'^ 

(A.l) 

■^Subject to the usual caveats of applying AdS/CFT techniques to QCD. 
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h equation: 



Ar^f {rB' + 4S - 4) (-6/ + + /) h' = 

C {Ar^CfhB' + IGr^Cfh {-Bf + C^+ /)) 
+4r^fhB' (/ - 3/^2) {-Bf + C^ + f) 

+e {-4r^fh''B' {-Bf + C^ + f)- Sr^C^h^ {-Bf + C^ + f)) 
{2r^Cf^hB' + 4r^C^f{r)h) 

+(t>'^ {r'^i-iB - l))f^hB' - 2r^iB - l)C^fh) 

+$'2 {-r^f^hB' - 2r^C^f%) - 2r^f^hB'^ 

+2r^//i,C'2 (C^ - 2{B - 1)/) - 8g2c(r)/i3^2^$ (C^ - 2{B - 1)/) 

-Aq^{B - (2(/? - 1)/ - C^) + 4g2//i^^'^' (C' - 2{B - 1)/) 

-4rV/iC'^ (-3(S - 1)CV + 2(S - Iff^ + C^) 

+24r2C(r)2/i (/ - h'') {-Bf + C^ + f) 

(A.2) 

B equation: 

4r'f{-Bf + C^ + f) B" = 

^r^CfB'C - 4r^h^ei3' {-Bf + + f) 
-Ar'^B' (/ + 3/1^) {-Bf + + /) + 2r'^CfB'(t)'<^' 
+$'2 (4^6^^ + c2 + /) - r^/^Z?') - 4r8/'S'2 

-r'^{B - l)fB'<t)'^ - Ar^{B - l)fC'^ 

(A.3) 

^ equation: 



Ar^f{-Bf + C^ + f) e = 

-4r^h^e^' {-Bf + C^ + f)- (/ + S/i^) + + /) 

-Sr^h''^ {-Bf + C^ + f)+ r\-{B - I)) f {r)^' <}>'^ + Aq\B - l)h^^<i>^ 
+2r^Cf^'<j>'^' - SqWm + V//i'<e*' - r Ve'*" 

(A.4) 



$ equation: 

4r^f{-Bf + C^ + f) = 

B' (4r2C/<^' - 4rV^^') + C {Ar^Cf^' - Ar^{B - l)f(j>') 
+8g2/i2^2^ + c2 + /) - 4r/i2^2^' + C2 + /) 

(4r (/ - 3/i2) {-Bf + C2 + /) - r(S - l)/,^'^) + 2rC/</)'$'2 
-r/2$'3 

(A.5) 
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C equation: 

Ar'^f {rB' + 4B - 4) {-Bf + + /) C" = 

e2</)$( - Sq^rh'^B' {{B-l)f+C^) + IGq^riB - l)Ch?C' 

+-i2q\B-l)h^ i-Bf + C^ + f)) 
+^^$2 {^q^rCfh^B' - 8q^r{B - \)fh^C') 

_ Ar'^h^B'C {-Bf + + /) - IQr^Ch'^B' {-Bf + + /) ) 
+C' (at^B' {{B - l)f (/ + 3/i2) + C2 (7/ - 3/i2)) 

_4r5/2B'2 _ g4^3(g _ + + /) ) 

+^'2 (8r^(B - 1)/C' + + /) - Sr^C/B' + + /)) 
+C'2 {Sr^CfB' - 16r^(B - 1)C/) 

+</.'2 (r4(-(S - 1))/B'C' - 4r^{B - l)CfB') 

+0'$' (2r^CfB'C' + Ar^fB' {-Bf + SC^ + f) 

+16r\B-l)f{-Bf + C^ + f)) 

(8gV(B - l)Ch^B' - 8q\{B - ifh^C) 
+ASr^CB' (/ - h^) {-Bf + + /) - 16r^Cf^B'^ - Ar^{B - l)fC'^ 

(A.6) 

(f) equation: 

Ar^f {rB' + 4B - 4) {-Bf + + /) (j)" = 

C'(b' {Ar^Cfcj}' - Ar^f^^') - lQr^{B - 1)C/.^' 
+16rV$' (_B/ + 2C2 + /)) 

{Sq^r^h^B' {-Bf + C^ + f) + i2q^r{B - l)h^ {-Bf + + /)) 

+f $ {l&q^{B - l)Ch^4)' - leg^C^/i^*') ) 

+B'(4r^.^' (3(B - 1)/ (/ + /i^) (/ - B/i^)) 

+r\-{B - l))f<f>'^ + {2r^Cf<p'^ - ISr^C/^) - rV'<^'*") 
+e {-Ar^h''B'cj>' {-Bf + C^ + f)- 16rW^' {-Bf + C^ + f)) 
+C'^ {Ar^Cf^' - Ar^{B - !)/</)') 

+$' {ASr^C (/ - /i^) + C^ + f)- 4r2(S - l)C/(/)'2) 

+^'2 {-8r\B - l)/0' ((S - 1)/ - C^) - 8r*C/$' {-Bf + C^+f)) 

-i2r^{B-l)f<P' {-Bf + C^ + f) 

+^2^2 ^g^2^^ _ i^c/i^^' - 8g2(e - lfh^(P') + Sr^C'^fcp'^'^ 
-Ar^Cf^<^'^ 

(A.7) 
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